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ABSTRACT 

The orbital decay and tidal disruption of a star cluster in a galaxy is studied in 
an analytical manner. Owing to dynamical friction, the star cluster spirals in toward 
the center of the galaxy. Simultaneously, the galactic tidal field strips stars from the 
outskirts of the star cluster. Under an assumption that the star cluster undergoes a 
self-similar evolution, we obtain the condition and timescale for the star cluster to reach 
the galaxy center before its disruption. The result is used to discuss the fate of so-called 
intermediate-mass black holes with > 10^ Mq found recently in young star clusters of 
starburst galaxies and also the mass function of globular clusters in galaxies. 

Subject headings: stellar dynamics — galaxies: star clusters — galaxies: starburst 


1. INTRODUCTION 


Tremaine, Ostriker, &: Spitzer (1975) studied a circular motion of a star cluster around the 
center of a galaxy (see also Binney &: Tremaine 1987, p. 427). The star cluster interacts with 
background stars of the galaxy via dynamical friction, loses its orbital kinetic energy, and thereby 
spirals in toward the galaxy center. Tremaine et al. (1975) approximated the mass density pb of 
the background stars with that of the singular isothermal sphere. 


PhiR) = 


AttGR^ ’ 


( 1 ) 


where R is the galactocentric radius, Ub is the circular velocity of the background stars, and G is 
the gravitational constant. The cluster mass M was set to be constant. The circular velocity of 
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the star cluster was set to be v\^. Then the rate of change of the galactocentric radius of the star 
cluster was obtained as 


= -0.4276 In Ab™, 
at Ub 


( 2 ) 


where InAb is the Coulomb logarithm. Solving equation (2) with the initial condition R = R\ 
t = 0, Tremaine et al. (1975) found that the star cluster reaches the galaxy center at the time 

t Aiecay; 


_ 1.169 i?fub 
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( 3 ) 


This timescale Tdecay has been often used to study the evolution of a star cluster. 


However, in practice, the star cluster is subject to the galactic tidal field. Stars belonging to 
the cluster have to lie within the tidal radius from the cluster center (Spitzer 1987, p. 101; see 
also Capriotti & Hawley 1996); 
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( 4 ) 


where My,{R) is the total mass of the background stars within the cluster’s galactocentric radius 
R, i.e., Mb(i?) = dvr p\yR^dR. If equation (1) is used for the mass distribution pb, equation (4) 
becomes 

Thus, with a decrease of the galactocentric radius R, the tidal radius rt becomes small. Stars are 
accordingly stripped from the outskirts of the cluster and become part of the background stars. 
Since the cluster mass becomes small, the dynamical friction becomes inefficient. In addition, if 
the star cluster is not sufficiently massive, it dissolves before it reaches the galaxy center. 


This effect of tidal stripping was incorporated in some past studies of the orbital decay of a 
star cluster, but they were numerical studies with limited ranges of the parameters. A more general 
analytical study is desirable. 

Here we conduct the analytical calculation for the first time. The star cluster is assumed to 
evolve in a self-similar manner (§2). Then the rate of change of the cluster mass is obtained as a 
function of the galactocentric radius of the star cluster (§3, eq. [14]). Using this mass loss rate 
together with equation (2) for the rate of change of the galactocentric radius of the star cluster, 
we derive the condition and timescale for the star cluster to reach the galaxy center before it 
dissolves (§4, eqs. [27] and [28]). We discuss the validity of our assumptions and the astrophysical 
applications of our result (§5). The latter discussion focusses on the fate of so-called intermediate- 
mass black holes with mass > 10^ Mq found in young star clusters of starburst galaxies and also 
on the mass function of globular clusters in galaxies. 
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2. SELF-SIMILAR EVOLUTION 


The description of a star cluster requires at least two independent internal parameters, e.g., 
the cluster mass M and the half-mass radius rh, in addition to external parameters such as the 
galactocentric radius R. To simplify our calculation, we consider a self-similar evolution, where 
one of the two internal parameters is exceptionally determined by the other. Here is an idealized 
self-similar model for a star cluster in a tidal field, which originates in the work of Spitzer (1987, 
p. 59; see also Capriotti & Hawley 1996). 

Recall that the crossing timescale Tcross — is much less than the other timescales 

relevant to the evolution of a star cluster. The star cluster is always nearly in virial equilibrium 
(Binney & Tremaine 1987, p. 211); 


E = K + W = -K = 


W 

T’ 


( 6 ) 


where E is the total energy, K is the kinetic energy, and W is the potential energy of the star 
cluster. Their definitions are 


K = 


M{v^) 

2 


and W 


M($) 

2 


(7) 


where v is the velocity of a star, is the gravitational potential at the position of a star, and (•) 
denotes an average over stars in the cluster. Spitzer (1969) observed that various spherical clusters 
approximately have the relation 



2GM 

5rh 


( 8 ) 


Then we have E = which yields 


dE _ GM"^ drh 2GM dM 
dt 5r^ dt 5rh dt 


This change in total energy dE is due to mass loss across the tidal radius rt and is thereby equal 
to the change in potential energy associated with the displacement of dM from the tidal radius to 
infinity. Thus we also have 

dE _ GM dM 
dt rt dt 

Equations (9) and (10) combine to yield 



din rh _ 2 ^2 _ dlnrt 
dt \ rt y dt 


( 11 ) 


where we have used M oc rj) (eq. [5]). If rh/rt = 1/3, this ratio remains the same. Since the 
tidal radius rt is determined by the cluster mass M and the galactocentric radius R (eq. [5]), the 
evolution of the star cluster is also determined by these two quantities. 
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3. MASS LOSS RATE 


The escape of the individual stars from the cluster across the tidal radius is due to two-body 
relaxation. This process is thereby represented by removing a fraction /esc of the stars every 
relaxation timescale Tj-eiax (Binney & Tremaine 1987, p. 490 and p. 523; Spitzer 1987, p. 52): 


dM 

dt 




"^relax 


-M. 


( 12 ) 


The half-mass relaxation timescale is used for Treiax (Spitzer & Hart 1971; Binney & Tremaine 1987, 
p. 514; Spitzer 1987, p. 40); 


_ 0.1386 M f rl 

Aeiax - J > 


(13) 


where In A is the Coulomb logarithm for the star cluster and m is the mean stellar mass. Us¬ 
ing equation (5), we rewrite equation (13) as Treiax = (0.05335/In A) x (r^/rt)^^^ x {MR/mv\^). 
Substituting this into equation (12), we obtain 


dM 

dt 


-18.741nA 


/ ?’h \ fesc-mvh 

\n) R 


(14) 


If rh/rt is constant, the mass loss rate in equation (14) becomes larger with a decrease of the 
galactocentric distance R. This is because the tidal radius rt, the half-mass radius and hence 
the relaxation timescale Ti-eiax become smaller. 


The escape probability /esc is approximated with the fraction of stars in a Maxwellian dis¬ 
tribution that have velocities exceeding the root-mean-square escape velocity (Uggc)^'^^ (Binney & 
Tremaine 1987, p. 490; Spitzer 1987, p. 51 and p. 57): 


47r 



(15) 


The escape velocity Vesc is defined at a position in the star cluster as 


vjcir) 

2 


+ <I)(r) 


GM 

rt 


(16) 


The average of equation (16) is taken over positions of stars in the cluster. 
(6)-(8), we obtain 

(vie) 

{v^) rt' 


Using also equations 


(17) 


Thus the escape probability is constant in a self-similar model. The ratio r^/rt = 1/3 leads to 
(Uggc)/(u^) = 7/3, which in turn leads to /esc = 0.07190. 
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ANALYTICAL SOLUTION 


With the initial condition R = Ri and M = Mi at t = 0, the behavior of a self-similar star 
cluster is studied by solving equations (2) and (14). The Coulomb logarithms InAb and In A, the 
radius ratio rh/rt, and the escape probability /esc are set to be constant. We introduce the friction 
timescale Tfric and the mass loss timescale rioss^ 

^ Rfvh 
4=0 ~ 0.4276 In AbGMi’ 


"Trie — 


R 


and 


—dR/dt 

M 


"noss — 


—dM/dt 


t=o 


18.74 In A/escm-Ub' 


(18) 

(19) 

rioss ’InAb \rt J GMf ’ 

The galactocentric radius of the cluster R, the cluster mass M, and the time t are normalized as 


We also introduce the ratio between Tfric and tiqss) 


Tfric In A f r\,\ 

a = JllS. = 43 . 83 - ' ' 


R M ~ ^ t ~ 

— = 72, — = M, and - =t. 

R'l M'l ’^fric 


Equations (2) and (14) are accordingly rewritten as 


dR _ M dM 

di R di 


a 

R' 


( 21 ) 


( 22 ) 


The initial condition is i?i = Mi = 1 at / = 0. 

Equation (22) yields MdM/dt = adR/dt. The solution of this equation is 

M^ = 2aR + l-2a. 

Using equation (23), we eliminate M from equation (22): 

R 


dt 

dR 


(2a.R + l-2a)V2' 


(23) 


(24) 


The solution of equation (24) gives the relation between the time t and the galactocentric radius 
of the star cluster R: 


t = 


(1 - 2a - a.R)(l - 2a + 2 aRyG - (1 - 3a) 
3a2 


(25) 


The relation between the time t and the cluster mass M is obtained by substituting equation (23) 
into equation (25): 

(3 - 6a - M2 )m - 2(1 - 3a) 


t = 


ba^ 


(26) 
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Figure 1 illustrates the evolutions of R and M as a function of t for a ^ 0 and a = 1/4, 1/2, 3/4, 
and 1. 


The condition for the star cluster to reach the galaxy center before its dissolution is M > 0 at 
i? = 0, for which equation (23) yields 


1 

a < - or 


"Hoss 2 


(27) 


The time at which the star cluster reaches the galaxy center is obtained by substituting R = 0 into 
equation (25); 


Tiecay 


Tdecay _ (1 - - (1 - 3a) 


Tfi-ic 3a^ 

If the parameter a is very small (a <C 1/2), equation (28) is reduced to 


(28) 


Tiecay 


1 a 

2 Tiecay 


T'ric "Trie \ 

2 V 3rioss/ ’ 


(29) 


In the limit a ^ 0, equations (28) and (29) reproduce the solution fdecay = 1/2 for no mass loss 
obtained by Tremaine et ah (1975, see our eq. [3]). The timescale fdecay for a = 1/2 is 2/3. Thus 
mass loss lengthens fdecay by a factor of 4/3 at most because significant mass loss occurs only at 
the latest stage of the orbital decay, i.e., in the vicinity of the galaxy center (Fig. 1). The cluster 
mass at the time when the star cluster reaches the galaxy center is obtained by substituting R = 0 
into equation (23): 

Mdecay = (1 - OI Mdecay = M (^1 - ' . (30) 

On the other hand, if the star cluster dissolves before it reaches the galaxy center, the corresponding 
time is obtained by substituting M = 0 into equation (26); 


_ 1 

^dissolve 

a 


1 

3a2 


or 


^dissolve Toss 



(31) 


The corresponding galactocentric radius is obtained by substituting M = 0 into equation (23): 


f^dissolve f n f^dissolve Ri I 1 r, 

2a V 2rfric 


(32) 


Therefore, the behavior of the star cluster is determined by the parameter a = ruic/rioss- If 
a <C 1/2, the mass loss is slow (2rfi.ic Toss)- The star cluster reaches the galaxy center before it 
loses much of its mass. If a 3> 1/2, the mass loss is fast (2rfric ^ Toss)- The star cluster dissolves 
before it moves significantly from its initial galactocentric radius. 
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5. DISCUSSION 
5.1. Model Assumptions 

We have assumed that a star cluster loses its mass via tidal stripping and two-body relaxation 
alone. This assumption is correct only if the orbit of the star cluster is circular and in the galactic 
plane. The star cluster otherwise suffers from gravitational shocks in passages through the galactic 
bulge or disk (Ostriker, Spitzer, & Chevalier 1972). These complicated cases are beyond the reach 
of our analytical calculation. 

We have assumed that the star cluster always has = 1/3. This ratio is unstable in our 

idealized model (eq. [11]). If ^h/rt > 1/3, the ratio increases still more. The star cluster dissolves 
at rh/rt = 4/5 where the escape velocity Vesc is zero (Capriotti & Hawley 1996; see our eq. [17]). 
If rji/rt <1/3, the ratio decreases toward zero. The ratio typical of globular and open clusters is 
rh/rt — 0.2 (Binney & Tremaine 1987, p. 26). Nevertheless, in practice, such star clusters do not 
undergo rh/rt ^ 0. The ratio is rather sustained to be nearly constant by heating due to binaries, 
a process that is not considered in our idealized model (Binney & Tremaine 1987, p. 543; Spitzer 
1987, p. 148). Since the mass loss rate for rh/rt = 0.2 obtained from equations (14), (15), and (17) 
is lower only by 12% than that for rh/rt = 1/3, our model with rh/rt = 1/3 is of practical use. 

There are other self-similar models for a star cluster in a tidal field. They are based on different 
assumptions. For example, H&on (1961) assumed the presence of an energy source at the cluster 
center and obtained rh/rt = 0.1446 and /esc = 0.045 (see also Spitzer 1987, p. 59).^ The energy 
source could represent the binary heating and sustains the small value of rh/rt as compared with 
that of our present model. It is possible to adapt these self-similar models into our calculation 
because the escape probability /esc is constant in any self-similar model (eqs. [15] and [17]). We 
nevertheless favor our present model, which is simple but still sufficient for our purpose. 

We have assumed that a star with velocity greater than the escape velocity immediately es¬ 
capes from the cluster. With somewhat improved treatments of the escape condition, the escape 
probabilities /esc close to ours were obtained by Fokker-Planck numerical simulations of Lee & 
Ostriker (1987) and Lee & Goodman (1995). 

The escape probability /esc is affected by two remarkable processes that are not considered 
in the above calculations. First, since stars with velocities slightly above the escape velocity can 
escape only through small holes near the Lagrangian points, they stay in the cluster for many 
crossing timescales. Some of the stars are scattered to lower velocities and become bound again. 
The escape probability is accordingly lowered (Baumgardt 2001). It might be better to consider 
that our /esc value tends to be an upper limit. Second, when the star cluster is young, mass loss 
in the course of evolution of the individual stars, i.e., stellar wind and supernova explosion, causes 


^Henon’s model was numerical. If we use our analytical approximation in equations (15) and (17), rh/rt = 0.1446 
yields /esc = 0.02006. 



shrink of the tidal radius and thereby enlarges the escape probability. This process is nevertheless 
unimportant in our study, which applies mainly to relatively old clusters. 

Taking account of the distribution of stellar masses as well as the evolution of stars and binaries, 
Portegies Zwart et al. (2002) conducted realistic A^-body numerical simulations of a star cluster in 
a static tidal field. While they did not vary the total number of stars in the cluster, they varied 
its initial density distribution and its distance from the galaxy center. The mass loss timescale 
Toss was found to scale with the relaxation timescale at the tidal radius (eq. [13] but with rt 
substituted for rh).^ Even if this is generally the case, our approach is still of use. We only have 
to set /esc = Teiax/Toss with r^/rt = 1. 


5.2. Astrophysical Applications 


The Coulomb logarithms In Ab and In A are of order 10^, and the circular velocity Ub is of order 
10^ km s“^ (Binney & Tremaine 1987, p. 423 and p. 427). Thus equations (18)-(20) are written 
as 


Trie = 5.316 X 10^ yr 


/In Ar 


V 10 


-1 


Vh 


100 km s“^/ \10® Mq 


Mi 


-1 


Ri 

1 kpc 


(33) 


= 1-387 X 10“ yr(!^) ‘ (^ 


3/2 


u. 


-1 


m 


1 Mp 


© 


-1 


Vh 


100 km s 


-1 


-1 


^0.07190 

Mi \ f R, 


10 ® Mq y V 1 kpc 


and 


„ = 0.03807 /. 

Toss VlnAb/ V 1/3 


m 


Vh 


1 Mq y V100 km s 


-1 


Mi 


10 ® Mq 


0.07190 
-2 


i?i 


1 kpc 


(34) 


(35) 


These formulae allow us to obtain the actual values of Tiecay) Afdecayj ©dissolve; and i?dissoive (eqs. 
[28], [30], [31], and [32]). The condition for the star cluster to reach the galaxy center before its 
dissolution, a < 1/2 (eq. [27]), is written as 


Mi > 2.759 X 10® Mq 


In A 


1/2 


/ ©h/©t \ 

InAb/ V 1/3 y 


m 


1 Mr 


0 


1/2 


©b 


100 km s 


-1 


Rj 

1 kpc 


-3/4 


1/2 


/esc \ 

0.07190 / 


1/2 


(36) 


^Within a static tidal field, the cluster mass decreases linearly with time (Spitzer 1987, p. 58). This behavior is 
reproduced by our mass loss rate (eq. [14]), if the galactocentric radius of the star cluster is set to be constant. 
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Figure 2 illustrates this condition as well as the values of rjecay and Tdissoive as a function of the 
initial mass Mi and the initial galactocentric radius Ri. 


5.2.1. Fate of Intermediate-Mass Blaek Hole 

Our result is used to discuss the fate of massive black holes (BHs) found in young star clusters 
of starburst galaxies. At the 2 ^m secondary peak of the starburst galaxy M82, i.e., an active site 
of star formation, there is a source of compact X-ray emission (Matsushita et al. 2000; Kaaret et 
al. 2001; Matsumoto et al. 2001). The observed strong variability implies that the source is an 
accreting BH. The observed luminosity of 10^^ ergs s“^ implies that the mass is greater than 10^ Mq 
if the emission is isotropic and its luminosity is below the Eddington limit. BHs of similar masses 
have been found also in other starburst galaxies, e.g., NGC 3628 (Strickland et al. 2001). They are 
called intermediate-mass BHs (Taniguchi et al. 2000) because their masses are in between those of 
stellar-mass BHs arising from stellar evolution and supermassive BHs found as the central engines 
of active galactic nuclei. The likely origin is successive mergers of massive stars and stellar-mass 
BHs that had segregated into the core of a star cluster (Taniguchi et al. 2000; Ebisuzaki et al. 2001; 
Portegies Zwart & McMillan 2002). Under a favorable condition, timescales of these processes are 
less than 10^ yr (Mouri & Taniguchi 2002b, c). 

Owing to dynamical friction, the host cluster of an intermediate-mass BH sinks to the galaxy 
center. Upon dissolution, the cluster releases the BH. If there are more than one such BH, they 
merge with each other. Thus the intermediate-mass BH could evolve to be supermassive and serve 
as the central engine of an active galactic nucleus (Matsushita et al. 2000; Ebisuzaki et al. 2001).^ 
If this is the case, intermediate-mass BHs are crucial to studying a possible evolutionary connection 
between a starburst and an active galactic nucleus. 

The host cluster is more massive than the BH and is affected by dynamical friction more 
significantly. Only if the host cluster reaches the galaxy center before it dissolves, the BH can 
reach the galaxy center in a reasonable timescale. If the host cluster dissolves before it moves 
significantly from its initial galactocentric radius, the timescale for the BH to reach the galaxy 
center is unreasonably large (Ebisuzaki et al. 2001). 

Young star clusters in starburst galaxies have masses up to 10® Mq (e.g., Mengel et al. 2002). 
Largest-mass clusters contain largest numbers of massive stars and thus are most advantageous 
to form intermediate-mass BHs. Such clusters are also most advantageous to bring the BHs to 
the galaxy center. Star forming regions of starburst galaxies generally have galactocentric radii of 


®Two remarks. First, a single intermediate-mass BH is sufficient to explain the lowest-luminosity active galactic 
nuclei, where the mass of the central engine is as low as 10^ Mq (e.g., NGC 4395; Filippenko & Sargent 1989). 
Second, even if there remains no available gas when the intermediate-mass BH reaches the center, the galaxy has 
a chance to emerge as that with active galactic nucleus in the future, since supply of a large amount of gas to the 
nuclear region is repetitive (e.g., Mouri & Taniguchi 2002a). 
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R ~ 10^-10^ pc (e.g., Telesco, Dressel, k, Wolstencroft 1993). If the mass is 10® M©, a star cluster 
there reaches the galaxy center within the timescale 10®-10® yr (Fig. 2). The timescale is 10® yr 
for the host cluster of the intermediate-mass BH in M82, which lies 170 pc from the galaxy center. 
These timescales are very small as compared with the lifetime of a galaxy, 10^® yr. 


5.2.2. Mass Function of Globular Clusters 


Our result is also used to discuss the mass function of globular clusters in galaxies. The 
globular clusters have a preferred mass scale. Their mass function has a sharp peak at (1-2) x 10® 
Mq (e.g., Harris 1991). A plausible explanation is that the mass function was initially wide but 
later modified by dynamical friction, tidal stripping, gravitational shocks, and so on (e.g.. Fall & 
Rees 1977; Okazaki k Tosa 1995; Fall k Zhang 2001). Although the shape of the initial mass 
function is uncertain, it is possible to demonstrate that the cluster mass 10® Mq is advantageous 
for a star cluster to survive from dynamical friction and tidal stripping (see also Lee k Ostriker 
1987; Capriotti k Hawley 1996). 

Suppose that star clusters were formed Tage yr ago over ranges of the mass Mi and the galacto- 
centric radius Ri. On our M\~R\ diagram (Fig. 2), the star clusters that have survived up to now 
are those above the line for the timescale Tage- This line is concave. If R\ is fixed, star clusters with 
too large masses have disappeared into the galaxy center due to dynamical friction. Star clusters 
with too small masses have dissolved due to tidal stripping. There is a initial mass Msurvive with 
which star clusters have survived for the widest range of R\. This mass corresponds to a ~ 1/2 
(solid line in Fig. 2), where the net effect of dynamical friction and tidal stripping is minimal. 
Eliminating Ri from equations (33) and (35) with Tage = Tdecay = '^dissolve = 2rfric/3 and a = 1/2, 
we obtain 


M«, 


\ / ^age 

100 km s“^/ \10^® yr 


/e: 


0.07190 


2/3 


(37) 


Thus, at the present age of globular clusters Tage — 10^® yr, those with Mj ~ 10® Mq are relatively 
abundant. Their present masses are also of order 10® Mq in most cases because mass loss is 
significant only at the latest stage of the orbital decay. 


This work has been supported in part by the Japanese Ministry of Education, Science, and 
Culture under grants 10044052 and 10304013. The authors are grateful to the referee for helpful 
comments. 
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FIGURE CAPTIONS 

Fig. 1.— Evolutions of R = R/R\ (a) and M = M/M\ (6) for a = ruic/rioss ^ 0 and a = 1/4, 
1/2, 3/4, and 1 as a function of t = t/rfnc. Here R is the galactocentric radius of the star cluster, 
R\ is its initial value, M is the cluster mass, M\ is its initial value, t is the time, Tf^ic is the friction 
timescale (eq. [18]), and rioss is the mass loss timescale (eq. [19]). The dotted lines show the R-t 
relation for the dissolution of the star cluster (a) and the M-t relation for the arrival of the star 
cluster at the galaxy center (6). 

Fig. 2.— Timescales Tdecay with which the star cluster reaches the galaxy center {dashed lines, eq. 
[28]) and Tdissoive with which the star cluster dissolves {dotted lines, eq. [31]). The abscissa is the 
initial value of the cluster mass M\ in units of Mq. The ordinate is the initial value of the cluster’s 
galactocentric radius R\ in units of pc. The other parameters are the same as in eqs. (33)-(36). 
The solid line corresponds to a = ruic/rioss = 1/2 (eqs. [35] and [36]) and divides the regions where 
the star cluster reaches the galaxy center before it dissolves {right) and the star cluster dissolves 
before it reaches the galaxy center {left). 















